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ABSTRACT 

Safety of the next-generation space flight vehicles requires 
development of an in-flight Failure Detection and 
Prognostic (FD&P) system. Development of such system 
is challenging task that involves analysis of many hard 
hitting engineering problems across the board. In this 
paper we report progress in the development of FD&P for 
the re-contact fault between upper stage nozzle and the 
inter- stage caused by the first stage and upper stage 
separation failure. A high-fidelity models and analytical 
estimations are applied to analyze the following sequence 
of events: (i) structural dynamics of the nozzle extension 
during the impact; (ii) structural stability of the deformed 
nozzle in the presence of the pressure and temperature 
loads induced by the hot gas flow during engine start up; 
and (iii) the fault induced thrust changes in the steady 
burning regime. The diagnostic is based on the 
measurements of the impact torque. The prognostic is 
based on the analysis of the correlation between the 
actuator signal and fault-induced changes in the nozzle 
structural stability and thrust/ 
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1. INTRODUCTION 

Stage separation failure may have various origins 
including e.g. failure of the accelerating/decelerating 
motors during the separation. The combined expected 
failure rate is relatively high and is the second most 
common cause of launch failures (Futron Corporation, 
2004). Despite of this fact a reliable in-flight diagnostic 
and prognostic system for this fault is currently 
unavailable. Below we report a progress in development 
of such system. 

We consider two-stage, vertically stacked space 
vehicle. The vehicle is propelled by the first stage booster 
at the first phase of the flight. After the rocket booster 
completes its mission, the stages separate and the second 
stage engine provides the propulsion during the second 
phase of the flight. Due to separation faults, the first stage 
can collide with the second stage engine nozzle and 
damage it. The hot gas flow may result in further 
degradation of the nozzle, reduction of the actual thrust, 
and significant side loads. 

Development of an in-flight FD&P system for this 
fault is a challenging engineering problem (Seife, 2003). 
The main difficulties stem from the fact that the 
phenomena underlying fault dynamics are highly non- 
linear in structural, thermal, and fluid dynamical domains. 
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Figure 2 Sketch of the separation failure geometry. Two 
fault induced separation trajectories are shown by black 
lines. Left path corresponds to sever impact: worst case 
scenario. Right path corresponds to a weak impact. 


Accordingly, the solution requires 3D analysis of the 
thermal/fluid/structure interaction in the supersonic flow. 
In addition the number of sensors available onboard is 
severely limited and the safe time window between the 
detectable onset of the fault and possible catastrophic 
failure is typically a few seconds. 

To overcome these difficulties we use a combination of 
physics based analysis and high-fidelity modeling. The 
analysis is performed for the worst case scenario and 
divided in three steps. In the first step (Sec. II) the 
structural dynamics of the nozzle during the separation 
failure is analyzed using finite element model and the 
theory of shells. It is shown that the nozzle buckles as a 
result of the impact. At the second step (sec. Ill) the 
pressure and temperature loads on the walls of the 
deformed nozzle in the hot gas flow are estimated 
analytically and numerically. It is shown that the most 
plausible consequence of the nozzle buckling is melting 
and disintegration of the thin nozzle walls in the region of 
the maximum plastic strain. High-fidelity simulations are 
employed to estimate the reduction of the actual thrust and 
side loads induced by the fault. In Sec. IV, we use the 
results of the preceding analysis to suggest the FD&P 
system of the stage separation fault based on the TVC 
actuator measurements of the fault-induced torque. 
Finally, we provide the discussion of the results and 
outline future direction of research. 


2. STRUCTUAL DYNAMICS 

A typical geometry of a separation failure is sketched in 
the Figure 1 . Fault induced by lateral velocity of the inter- 
stage results in a separation trajectory that crosses the 
nozzle. The greater is the lateral velocity the larger is the 
overlap of the inter- stage and he nozzle and more severe is 
the consequence of the impact. The goal of the FD&P 
system is to detect the impact, estimate its strength, and to 



Figure 1 Sketch of the nozzle as a conical shell defined by 
the coordinates s and 0. 

predict the fault-induced torque of the damaged nozzle 
after engine ignition. 

is to detect the impact, estimate its strength, and to predict 
the fault-induced reduction of the actual thrust and side 
loads. 

For the typical parameters of the nozzle extension of 
the liquid motor the nozzle dynamics under the interstage 
impact can be considered based on shell theory [1]. A 
simplified geometry of the nozzle extension 
corresponding to this model is shown in the Fig. 2. We 
consider dynamics of the second stage engine nozzle 
extension using Donnell's equation [1] 


dw 


w 
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dt 
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where w transverse displacement, N- an Airy type of 
stress function, D and G are the bending and shear 
stiffness c - damping coefficient, h - thickness of the 
conical shell, R is the mean nozzle radius, f(t,s,0) local 
external pressure acting on the nozzle from insterstage 
and 
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are operators in cylindrical coordinate system, 6 is the 
angular coordinate. 

We consider boundary conditions for the trancated 
conical shell, short end of which is clamped and free large 
end is free. To find the eigenfunctions and 
eigenfrequencies of the problem we consider the next 
eigenvalue prolem 
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where 
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The eigenfunctions w mn (s), N mn (s ) satisfying Eq. (2), and 
natural frequencies of the cone can be found out by 
Galerkin's or Rayleigh— Ritz methods. Expanding the 
radial displacement w and Airy function N in the series of 
the shallow shell eignemodes 

w(t,s,0) = cosine), 

mn 

N{t,s,0) = Yd mn (t)N mn {s) cosinG) 


where n is the number of waves in the circumferential 
direction and using the property that normal modes with 
different indexes are orthogonal we obtain equations 

governing the amplitudes C mn as 


dV dr 

2HL + 2j^HL + (tf> 2 ) C = f ( t ) (3) 
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Amplitutes d mn differe form c mn solution by the constant. 
Solution of the equation (3) with the intial condition 

can be presented by Duhamel’s integral [2] 


J / mn (r)sin« mn (r-r)Jr 



This approach makes it possible to calculate elastic stress- 
strain distribution in whole nozzle during the impact. For 
the impact time of the order of 50 msec we expect that the 
frequency range of the eigenmodes dominating the 
response dynamics will be limited by 10 Hz. 

The linear elastic approach was justified by estimation 
of buckling stresses during impact. The external load 
f(t,s,6) on the local surface can lead to sufficiently large 
stresses and nonlinear dynamics determined by plasticity. 


The external buckling stresses corresponding to the lateral 
pressure of the conical shell was estimated by formula [1] 


Vforg h_(h V /2 
a ~ 9 ^ v 2f l{r) 



( 4 ) 


where E is the Young’s modulus of the cone, L is the 
length of the cone which is of the same order with radius 
(L~R). Considering £=2xlO u Pa , h/R~ 2xl0' 3 we have the 
buckling pressure 2xl0 7 Pa , which is of the same order, 
observed in the nozzle during impact simulations. This is 
why Finite Element model was developed. 

To analyze response of the nozzle to the impact we 
build a finite element model of the nozzle extension in 
ABAQUS. Because the nozzle shell is much thinner then 
the impacting part the later is modeled as a rigid body. 
We fix the nozzle extension at its base and use a dynamic 
explicit mode with general contact properties to simulate 
the impact. Typical results of the simulations are shown in 
the Figure 3. It can be seen from the figure how the 
nozzle is damaged during the impact. One of the main 
conclusions of the simulations is that the nozzle buckles 
as a result of the impact. The buckling stresses found in 
the numerical simulations confirm the values obtained 
with analytical estimations above. 

Next we estimate the impact-induced torque on the 
nozzle. The results are shown in the Figure 4 for three 
different values of the overlap of the nozzle extension 
with the impacting cylindrical shell at the initial moment 
of the impact. 

It can be seen from figures that the strength of the 
impact can be inferred from the analysis of the torque 
signal and specifically from the amplitude of the first 
impulse. In addition the lateral relative velocity of the 
separation is usually associated with the malfunction of 
the accelerating/decelerating motors. As a result axial 
velocity is also affected by the fault and there is a fault- 



Figure 3 Results of the simulations of the impact. The 
bended and buckled areas of the nozzle are shown by the 
arrows. 
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Figure 4. Impact-induced torque on the nozzle for three 
different overlaps of the nozzle extension with the 
impacting cylindrical shell at the initial time of the 
impact. 



Axial Position 


Figure 5. Radial pressure distribution in quasi-two- 
dimensional approximation is shown for four different 
locations along the nozzle axis. Radial pressure 
distribution (circles) as compared with the results of 
high-fidelity simulations (dashed line). 


induced delay time between the separation initiation and 
the impact instant. The more sever is the fault the larger is 
the delay. Therefore correlation between the amplitude of 
the impact torque and the impact delay time can serve as 
an additional characteristic dynamical signature of the 
separation fault that helps to infer the fault induced 
damage. 

Our analysis shows that the damage can be 
characterized by two key factors: (i) the area S of the dent 
and (ii) the angle of the bend. For the separation failure 
geometry shown in Figure 2 the area of the damage can be 
as large as 0.5 m 2 and the angle of the bending of the 
nozzle vary from 30 to 50 degrees. It will be shown in the 
following discussion that the structural stability of the 
nozzle and fault induced changes of the thrust are mainly 
determined by the area and angle of the bending. 


consider a one dimensional isentropic model of the nozzle 
flow for given transient time traces of the pressure p 0 {t ), 
mass flow rate J(t), and temperature T 0 (t) in the main 
combustion chamber. The exit pressure p ex (t ) and 
temperature T ex (t) can be written in the form 


-T 

-*■ n 
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We note that the Mach number at the exit M ex = u ex /c 0 
does not depend on the flow variations and is determined 
by the nozzle geometry and the specific heat ratio y~ cjc v 
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3. TRANSIENT LOADS ON THE NOZZLE 

To estimate further degradation of the damaged nozzle in 
the hot gas flow one has to determine pressure and 
temperature loads on the nozzle walls caused by the flow. 
To do so we use combination of high-fidelity simulations 
and analytical estimations. Analytical estimations are 
based on the analysis of the flow parameters at the nozzle 
surface in the presence of the dent. To determine transient 
parameters of the flow a quasi-two-dimensional model of 
the flow introduced. This model is verified by comparison 
with static high-fidelity simulations. The transient loads on 
the wall are then introduced into the finite element model 
of the nozzle to analyze coupled thermal-structural 
dynamics of the damaged nozzle in the hot gas flow. 

3.1 Quasi-two-dimensional approximation 

To estimate the pressure and temperature loads on the 
damaged nozzle a quasi-two-dimensional model of the 
transient flow is introduced as follows. Let us first 


where S t is the nozzle throat area and S ex is the nozzle exit 
area. Thus time variation of the thrust is completely 
determined by the time variation of the chamber pressure 
according to the following expression 

F lh (t)=^p (] (t)S,M a . (7) 

To determine pressure and temperature loads on the 
nozzle wall we use the results of static high-fidelity 
simulations. The radial distribution of pressure obtained 
in high-fidelity simulations is shown in the insert of 
Figure 5 by open circles. We next require that: (i) for 
every location along the nozzle axis the mean value of the 
radial pressure distribution coincides with the values 
obtained using one-dimensional model (5); (ii) the wall 
pressure coincides with the results of the high-fidelity 
simulations; and (iii) the values of the vacuum thrust 
calculated using Eq. (7) correspond to the vacuum thrust 
obtained by integrating pressure on the nozzle walls. 
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Figure 6. Oblique shock wave with deflection angle a 
near the bending corner of the damaged area is shown by 
the dashed bended lines. Parameters of the incident and 
aftershock flaws are indicated by sub-indexes 1 and 2 
respectively. 

These conditions can be satisfied using polynomial fit of 
the high-fidelity radial pressure distribution and requiring 
that mean value of this polynomial approximation is given 
by Eq. (5). A quasi-two dimensional distribution of 
pressure in the nozzle obtained in such approximation is 
shown in the Figure 5. Similarly we obtain quasi-two- 
dimensional distribution of the flow temperature. The 
advantage of such approximation is that it allows us to 
obtain accurate representation of the transient nozzle flow 
parameters near the nozzle wall at the location of the 
damage. 

Next, we apply this approximation to analyze the shock 
structure near the damaged wall and to estimate the 
pressure and temperature loads on the damaged wall 
induced by the aftershock flow. 


3.2 Temperature and pressure loads 

To estimate temperature and pressure loads on the wall we 
recall that whenever a supersonic flow is “turned into 
itself’ as shown in Figure 6 an oblique shock wave will 
occur (Shapiro, 1953; Andreson, 2001). The analysis 
shows that for deflection angles a as large as 50° and 
relevant gas parameters the flow at the corner will 
correspond to the detached curved shock wave as shown in 
Figure. Accordingly, estimation from above for the 
pressure and temperature loads on the damaged wall can 
be performed using equations for the normal shock wave 
in the form 


i + Rr + i)/2K 

2 yM?-(y-l)/2 ’ 


= 1 + 
P i 


2y 

(r+ 1 ) 


K 2 - 1), 


r 2 i+[(r-i)/2]M, 2 

Ti l + [(r-l)/2]M 2 ' 


( 8 ) 



Figure 7. Distribution of the velocity of the gas flow in 
the damaged nozzle. The damaged area and the 
formation of the detached oblique shock wave can be 
seen in the left upper corner of the nozzle. 

To estimate temperature and pressure loads from 
below we use the smallest deflection angle a = 30° of the 
damaged area obtained in structural dynamics 
simulations, substitute in the first Eq. (8) the values of the 
Mach numbers M 7 and M 2 with the values of Mach 
numbers normal to the shock, and find the shock wave 
angle by solving the equation 


tan a- 2cot /? 


M 2 sin 2 {3-1 
M 1 2 (^+cos2^) + 2' 


( 9 ) 


The resulting pressure loads on the damaged wall 
range from approximately 1 atm (weak oblique shock 
with deflection angle a = 30°) to several atmospheres 
(normal detached shock). The corresponding temperature 
on the wall varies from 2200 K to 3700 K, which exceeds 
the melting temperature T md ~ 2150K. 

This analytical estimation were confirmed by the high- 
fidelity simulations of the gas flow in the damaged nozzle 
using FLUENT. An example of the velocity distribution 
obtained in these simulations is shown in Figure 7. The 
formation of the detached oblique shock wave can be 
clearly seen near the dent corner of the nozzle. 

Note that the obtained pressure and temperature loads 
may result in yielding and melting away of the nozzle 
dent as will be discussed in the following section. We 
now estimate the time scale required for yielding, melting, 
and burning through the damaged area of the nozzle. 


3.3 Fault-induced transient torque 

It was shown in the previous section that the temperature 
loads on the damaged wall exceed the melting 
temperature. Accordingly, the transient dynamics of the 
effect of the flow on the damaged wall will be mainly 
determined by the time required to heat the wall to the 
melting temperature. Simultaneously with nozzle heating 
the rise of the pressure load during the engine ignition 
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Figure 8. Dynamics of the temperature of the internal 
(blue solid line) and external (black dashed line) surfaces 
of the damaged nozzle. Red horizontal line shows 
melting temperature on the nozzle material. 


transient will exceed yielding limit required for reversed 
buckling just before the nozzle will begin to melt. 

To estimate the characteristic time scale of the 
corresponding dynamics of heating and buckling of the 
damaged nozzle, we use time traces of the ignition 
transient for the temperature and pressure reported by 
(Wang and Guidos, 2009). According to their results the 
pressure and temperature jump of the nozzle flow occurs 
1.7 sec after ignition when the nozzle is uniformly filled 
with the supersonic flow. The subsequent rise of the flow 
parameters to the nominal steady values occurs on the time 
scales from 0.5 to 1.0 sec with the fastest rising time (~ 
0.5 sec) corresponding to the temperature change. 

We use ABAQUS to integrate equation of structural 
dynamics of the damaged nozzle under transient pressure 
load coupled with the solution of the heat equation for the 



Figure 9. Torque induced by the transient dynamics of the 
aftershock temperature and pressure loads on the 
damaged nozzle area S - 0.15 m 2 . The inset shows the 
change in the shape of the damaged area 0.18 sec after the 
pressure jump. The initial shape is shown by the gray 
contour; the modified shape is shown in colored contour. 


wall temperature T where boundary conditions correspond 
to the heat flux from the gas flow to the internal surface of 
the nozzle with the heat transfer coefficient 77 and zero 
flux on the external surface of the nozzle. 

The temperature dynamics are shown in Figure 8. It 
can be seen from the figure that the temperature of the 
external side of the thin nozzle walls approaches melting 
temperature T md in approximately 0.23 sec. 

To estimate the characteristic time of the heating of 
the nozzle wall analytically we use Bartz' approximation 
(Bartz, 1965; Incropera and DeWitt, 2002) for the heat 
transfer from the gas flow to the nozzle wall. Since the 
nozzle wall is very thin the analysis is performed in ID 
using the following equation for the dynamics of the 
nozzle temperature T(x,t) 


r = 


K 


^ 'met P nu 

dT(x,t ) 


-K 


-K 


dx 

dT(x,t ) 


dx 


x 1 >x>x 0 , t>t 0 


\X-Xq 


= 0, T(x,t 0 ) = T 0 , 


( 10 ) 


Here internal surface is situated at x = x 0 and external 
surface is situated at v = xj. The solution of this equation 
gives an estimate for the characteristic melting time 
approximately 0.2 sec in good agreement with numerical 
results. 

However, even before the melting begins the pressure 
load on the wall exceeds the critical value for reversed 
buckling. Indeed, the critical stress can be estimated as o 
- pS/h 2 where p is the aftershock pressure, S is the 
damaged area, and h is the nozzle wall thickness. For h - 
5 mm, 5-0.5 m 2 , and given pressure transient time-traces 
the critical pressure p - 7 x 10 4 Pa is reached in 
approximately 0.16 sec after the pressure jump (i.e. 
approximately 1.85 sec after engine ignition). At this 
moment the reversed buckling of the damaged area of the 
nozzle is initiated by the pressure loads induced by the 
nozzle flow. The damaged nozzle wall will begin to melt 
down approximately 40 msec later. So in approximately 2 
sec after the engine ignition the damaged area of the 
nozzle will burned through by the gas flow. 

To analyze this sequence of events in more details we 
solve numerically the coupled temperature-displacement 
problem in ABAQUS for the relevant parameters of the 
flow (Hill and Peterson, 1992). The results of these 
simulations are shown in Figure 9. The analysis reveals 
that the reversed buckling begins at approximately 0.16 
sec after the pressure jump and the dent is completely 
buckled back 0.25 sec after the pressure jump. 

Thus our analysis shows that under plume loading the 
dent causes a reaction torque and temperature load on the 
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Figure 10. (top) The reduction of the actual thrust in the 
notched nozzle as function of the notch area S n . Different 
colors correspond to different meshes of the model. The 
eye guiding lines approximate the mean value (dashed 
lines) and the error of estimation (solid lines), (bottom) 
The side load developed in the damaged nozzle as a 
function of the notch area. The results of the numerical 
calculations are shown by magenta rhombus, the green 
and yellow circles show results of analytical estimations. 

nozzle until the damaged area yields and melts away, 
leaving a notch in the nozzle wall. 

4. QUASI-STEADY THRUST 

Finally we estimate reduction of the actual thrust and the 
side load in the notched nozzle in quasi-steady regime. 
The results of the numerical simulations for different 
values of the S n are shown Figure 10 (top). The thrust was 
calculated as a product of the velocity and the mass flow 
rate the at the nozzle exit surface. The analysis shows that 
the reduction of the actual thrust and the side load caused 
by the notch are both linearly proportional to the notch 
area S n . 

The results of the numerical calculation of the fault- 
induced torque are shown in Figure 10 (bottom), together 
with the results of analytical estimations. The later 
estimations were obtained in the point force 
approximation. It was assumed the force on the wall in 
radial direction is compensated everywhere except at the 
notch area and symmetrically situated undamaged area of 
the nozzle S n0 . The force on the notch was equal to zero 
and the force on the undamaged area was calculated as a 
point force and equal to the integral of the nominal 
pressure on the nozzle wall over S n0 . The estimations 



Figure 11. Transient torque induced by impact and the 
flow. 

shown by the green circler correspond to the assumption 
that there is no asymmetry in the axial load. The 
estimations shown by the yellow circles take into account 
the asymmetry of the load in axial direction. The results 
of the analytical estimations are in a good agreement with 
the results of numerical simulations. 

5. CONCLUSIONS 

The analysis of the torque induced by the impact and by 
the ignition transient loads on the damaged nozzle 
reported above allows one to build the FD&P system for 
in-flight detection of the stage separation failure and 
prediction of the fault induced thrust. 

As it was shown in Sec. II the impact induced torque 
provides a convenient tool for detection of separation 
failure and inferring the size of damage via measurements 
of the amplitude of the first peak of the impact torque. In 
addition, correlation between the amplitude of the torque 
and the time delay between the separation ignition and the 
impact signal can be used to improve the accuracy of the 
detection algorithm. It was shown further in Sec. Ill that 
fault induced damage of the nozzle will result in 
significant ignition transient torque imposed by the flow 
on the buckled area of the nozzle surface. The value and 
sign of this torque are strongly correlated with the 
strength of the impact. In particular, flow will induce 
torque in the same plane of separation in the opposite 
direction compared to the impact induced torque. It was 
also shown in Sec. Ill that the heat and pressure loads 
generated by the oblique shock wave near the dent corner 
of the damaged nozzle will result in reverse buckling of 
the damaged area with its subsequent melting and 
disintegration. The time scale of this process was 
estimated to be approximately 0.23 sec. As a result the 
torque induced by the flow will change its sign one more 
time and will approach quasi-equilibrium value 
proportional to the damaged area. For convenience we 
have summarized the torque time-traces at different stages 
of the separation failure in Figure 11. The in-flight 
analysis of these time-traces can be conveniently 
performed within a Bayesian inferential framework 
similar to the one introduced in our earlier work 
(Luchinsky, 2009; Osipov, 2007) for the in-flight 
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diagnostic and prognostic of the case breach fault in Solid 
rocket motor. The corresponding work is currently in 
progress. 


NOMENCLATURE 


gas density 
gas pressure 
gas temperature 
gas velocity 
sound velocity 
thickness of the nozzle wall 
Mach number, M = u/c, M 0 = u/c 0 
specific heat for the constant volume 
specific heat for the constant pressure 
ration of specific heats; y= C p /C v 
nozzle radius 

cross-sectional area of the nozzle 
area of the damaged part of the nozzle 
thrust produced by gas flow 
side thrust 

mass flow density [ j = up] 
heat transfer coefficient 
heat flow from the gas to a hole wall 
melting temperature point 
temperature of metal surface burning 
specific heat of case metal 
density of case metal 
thermal conductivity 
dynamical viscosity of hot gas 
vertical deflection of the nozzle surface 
bending stiffness 
shear stiffness 
nozzle radius 
nozzle wall thickness 
nozzle surface coordinates 
Airy type of stress function 
deflection angle 
shock wave angle 

Subscripts: 
in initial states 

ex parameters in the nozzle exit 

t parameters in the nozzle throat 

0 parameters in the chamber 


P 

P 

T 

u 

c 
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Cy 

Cp 
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F th 

F s 

j 

0 
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T me l 

T b 

Fmet 

Pmet 
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